Driving force for turbulent transition: Dynamical instability or statistical preference? 



Sergei F. Chekmarev 1 '^ 

1 Institute of Thermophysics, SB RAS, 630090 Novosibirsk, Russia and 
2 Department of Physics, Novosibirsk State University, 630090 Novosibirsk, Russia 

(Dated: January 13, 2013) 

A simple statistical model is proposed, which considers the turbulent flow as a collection of 
localized structures (eddies) . The model successfully predicts the transition to the turbulent regime 
at high Reynolds numbers and shows agreement with other known properties of turbulent flows. 
Comparison of the assumptions underlying the model and Navier-Stokes equation suggests that 
the driving force for the turbulent transition is basically the tendency of the system to occupy a 
statistically dominant state; the instability of the flow plays a role of the mechanism for finding this 
state. 
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The transition from laminar to turbulent fluid flow, 
originally studied by Reynolds in his seminal work [lj, 
is still in focus of many researchers. It is a very com- 
plex phenomenon, which has a variety of scenarios and 
is generally dependent on the flow environment 0-13. In 
the simplest cases, the transition is related to the dy- 
namical instability of a laminar flow at high Reynolds 
numbers, and this concept is well supported by the anal- 
ysis of the solutions of hydrodynamic equations [4, 5]. 
The transition leads to the formation of localized spatial 
structures (eddies) filling the flow field, which are consid- 
ered to be inherent elements of turbulent flows [6|-|9(. It 
is thus of interest to look at this phenomenon from an- 
other, "pure" statistical viewpoint, i.e. to see how much, 
if so, the structured turbulent flow is statistically more 
favorable than the unstructured laminar flow. For this, 
a simple model is proposed, which considers the turbu- 
lent flow as a collection of localized structures of vari- 
ous size and compares the statistical weights of different 
collections. Despite its simplicity, the model leads to 
meaningful results: it successfully predicts the transition 
to the turbulent regime at high Reynolds numbers and 
shows reasonable agreement with other known properties 
of turbulent flows. The model is based essentially on two 
assumptions, i.e. that the particles (atoms or molecules) 
constituting the system are identical (indistinguishable), 
and some elementary volumes of the phase space exist in 
which the states of the particles are uncertain. Since the 
same assumptions are involved in the derivation of hydro- 
dynamic (Navier-Stokes) equations, the obtained results 
suggest that the onset of turbulence is basically due not 
to the instability of the laminar flow but the tendency 
of the system to occupy a statistically dominant state; 
the role of the instability is then reduced to performing 
a structural rearrangement of the flow to find this state. 

Let us consider a system of N identical particles. Sup- 
pose that the particles can form localized structures of iVj 
particles, with the number of such structures being M». 
These structures are supposed to be similar to eddies in 
that the particles execute a concerted motion. Further, 
assume that each i structure consists of N%/n elemen- 
tary cells (each of n particles), such that the behavior 



of the particles in the cell is uncorrelated and random. 
Then, assuming that the interaction between the local- 
ized structures is small, the total number of states can 
be estimated similar to that for an ideal-gas system 



r = 



EL [(n^/^Ni/n)^ MA 

Assuming n, Ni/n and Mi to be large enough, so that 
the Stirling approximation x\ ks (x/e) x is applicable, and 
taking into account that ^\ NiMi — N, this equation 
rewrites as 



{N/n) N e N / n 



(1) 



With Eq.([T]), it is possible to calculate the most 
probable size distribution of the structures, Mj = 
Mi(Ni), which maximizes the number of states T and, 
correspondingly, the entropy S = lnT. This dis- 
tribution reduces the phase space of the system to 
a subspace of much lower dimension similar to how 
the phase space shrinks when the system reaches the 
(strange) attractor in the transition to turbulence 0, 
111 ]. Applying the Lagrange multiplier method to 
the entropy functional, i.e. varying the expression 
[-(MiNi/n) HNjn) - MMM l /e)}+a'J2 l M l N t + 
j3' Yli MiEi with respect to Mi at two conservation con- 
ditions 



and 



J2 MiNi = N 



MiEi = E 



(2) 



(3) 



where Ei is the kinetic energy of i structure, and E is 
the total kinetic energy, one obtains 



Mi = 



N 



^a'Ni+P'Ei 



(4) 



where and a' and /3' are the Lagrange multipliers. 
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The dependence of Ei on Ni is specific rather than 
universal; it is clearly different, e.g., for homogeneous and 
shear-layer turbulence. To be definite, suppose that the 
kinetic energy of fluctuations per unit mass increases with 
the distance r as in the Kolmogorov theory of turbulence 
for the inertial range of scales [7H9|, i.e. 



e(r) 



r 2/3 



(5) 



Then, the kinetic energy of the structure of Ni particles is 
Ei ~ Jq* e{r)r 2 dr ~ r] 1 ^ ~ N^ 1 ^ 9 . Also, it is convenient 
to pass from Ni to the number of the elementary cells in 
the localized structure qi — Ni/n. With these changes, 
Eq. becomes 



Mi 



(6) 



where a and /3 are new constants that are determined 
from the equations of conservation of the total number 
of particles [Eq. ((2J)] 



N/n = M i<li 
i 

and the kinetic energy [Eq. (J2J)] 

E/n n / s = J2M iq l 1/9 



(7) 



(8) 



Given the values of n, N and E, the constants a and 
/3 can be calculated by generating a set of qi and solving 
Eqs. and © with M l from Eq. ©. Another possi- 
bility is to vary a and j3 to obtain N/n and E/n 11 ^ 9 as 
functions of a and /3. 

The parameter N/n can be associated with the 
Reynolds number ReL = VL/v, where V and L are, 
respectively, the velocity and linear scale characterizing 
the system as a whole, and v is the kinematic viscosity 
of the fluid. Since the particles are identical, the 6N- 
dimensional phase space is reduced to the 6-dimensional 
single-particle space. Then, the linear size of the ele- 
mentary cell, in which the behavior of the particles is 
uncorrelated, can be taken as the product of the char- 
acteristic velocity and length at which the state of the 
particle is uncertain. For a gas fluid, it is the product 
of the molecular thermal velocity c and the mean free 
path A (l2j |. Then, the phase volume of the cell is (cA) 3 , 
with the number of particles in the cell being propor- 
tional to this volume, n ~ (cA) 3 (the density of the fluid 
is assumed to be constant). Correspondingly, the total 
number of particles TV can be considered to be propor- 
tional to the total phase volume for the system (VL) 3 , 
i.e. N ~ (VL) 3 . Taking into account the kinetic theory 
expression for the gas viscosity v ~ cA [l3j |. one obtains 



N/n = a(VL/uf 



aRe? 



(9) 



where a ~ 1. For a liquid, the linear size of the elemen- 
tary cell can be taken to be c 2 r, where c 2 presents the 




FIG. 1: Average number of the elementary cells in the local- 
ized structures (q) as a function of Re^; a and ft vary from 
-1.5 to 1.5. The up and down triangles correspond, respec- 
tively, to 1 < qi < 50 and 1 < qi < 100. 
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FIG. 2: The values of the parameters a and /3 at the left 
and right boundaries of the transition region of Fig. [T] for 
1 < qi < 100; a and ft vary from -1.5 to 1.5. 



fluctuations of the (molecular) kinetic energy per unit 
mass, and r is the mean residence time of the molecule 
in a settled state. According to the elementary kinetic 



theory [14] . the viscosity of liquid is v ~ err. We thus 
arrive at Eq. © again. In what follows, we assume for 
simplicity that a — 1 in Eq. i.e. 



Re L = (N/n) 



1/3 



(10) 



Figure [T] shows the average number of the elementary 
cells in the localized structures (q) = J2i Qi^i/^2i as 
a function of Re^ . The parameters a and (3 vary indepen- 
dently from a mm = /3 min = -1.5 to a max = ^ max = 1.5 
(every point thus represents a combination of a and /3), 
and 1 < q. t < g max (g max = 50 and 100). As g max in- 
creases, the transition region slightly extends to larger 
values of Re^, but the distribution of (q) at small to 
moderate values of Re^ practically does not change. It is 
seen that the size of the structures drastically increases 
from (ij) « 1 at Rej, ~ 1 to (q) w g max at ReL ~ 10 5 . 
Also, the dependence of (q) on Rez, is not unique, i.e. the 
same values of (q) are observed in a broad range of Re/, . 
At the left boundary of the transition region /3 = (3 max , 
whereas a increases with (g), and at the right boundary 
ol = a max , and (3 increases with (q) (Fig. In the 

given case the negative values of a and (3 are ineffective 
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FIG. 3: Average number of the elementary cells in the local- 
ized structures (q) as a function of Rez,; 1 < qt < 100. The 
solid and empty triangles are for a and /3 varied, respectively, 
from -1.5 to 1.5 and from -1.75 to 1.75. 
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FIG. 4: Characteristic linear size distributions of the localized 
structures. The dashed line curves show the corresponding 
Gaussian distributions; 1 < qt < 100, and a and /3 vary from 
-1.5 to 1.5. 



(because the pre-exponential factor in Eq. © decreases 
with qi very fast), however, as the interval of the varia- 
tion of a and ft increases, a becomes negative along all 
the left boundary. Figure [3] illustrates how the distribu- 
tion of Fig. Q]at qmax = 100 changes with the extension 
of the range of variation of a and (3. It is seen that the 
left boundary shifts to the lower values of Rei, and the 
right boundary to higher values of Rei,. The considera- 
tion of a wider range of variation of a and /3 shows that 
the Reynolds numbers characterizing the right boundary 
increase with a max as Rei, ~ exp(Aa max (q)) (A w 0.06), 
and those of the left boundary decrease with /3 max as 
Rei ~ exp(£?/3 max (<7)) (B w —0.13). The latter, however, 
continues until Re_i ~ 10 (/3 max ~ 2.1) is reached; after 
that the left boundary practically does not move with the 
increase of /3 max . Thus, the parameters /3 max and a max 
determine, respectively, the left and right boundaries of 
the transition region. 

Figure Q] suggests that the transition region from lam- 
inar ((q) ~ 1) to turbulent ((q) 3> 1) flow is between 
Rei, ~ 10 2 and Rei ~ 10 5 , which corresponds to the 
experimental and simulation results in hydrodynamic 
turbulence surprisingly well [H-13]- Moreover, the finite 
width of the transition range at a fixed value of (q) indi- 
cates that the transition can occur at different Reynolds 
numbers, which is also characteristic of turbulent flows. 
In his original work [l| , Reynolds has found that the criti- 
cal Re number for the pipe flow varied from Rei, w2x 10 3 
to Re^ ~ 1.3x 10 4 depending on the inlet conditions. The 
lower bound of the stability of the laminar flows, which is 
observed at large disturbances of the inlet flow, has been 
confirmed in many experimental and theoretical works 
(see, e.g., Refs. [l5| and [lfjj], respectively, and papers 
cited therein), and the upper bound, which is achieved 
in carefully controlled conditions, has been found as high 
as Rei w 1 x 10 5 [17j. In the present model, the above 
mentioned lower bound of the stability corresponds to 
/?max ~ 1-5, and the upper bound to a max « 1.75. 

Similarly to (N/n) 1 / 3 [Eq. ([TO])], the quantity qj /3 = 



(iVj/n) 1 / 3 can be considered to be the Reynolds number 

1 /3 

of i structure, i.e. Re^ — q i . According to the Kol- 
mogorov theory , Re^ ~ 1 characterizes the dissipation 
scale 77. Let us rewrite q t — N/n as (Ni/N) x (N/n) 
and take into account that according to Eq. (JS|) Vi/V ~ 
(k/L) 1 / 3 , so that Ni/N - [vtk/iVL)} 3 ~ (k/L)\ Then, 
the equality qi ~ 1 gives the well known result rj/L ~ 
Re^ 3 ^ 4 (see, e.g., Q). Since in the Kolmogorov theory 
Re^ = r/Vrf/v ~ 1 (v n is the velocity increment at the 
dissipation scale), it also follows that the linear size of 
the elementary cell could be determined as rjv^. 

Figure [4] shows the probability density distribu- 
tions of the linear sizes of the structures Mi(li/rf) = 
Mi(qi)dqi/d(li/rj) for three characteristic values of the 
Reynolds number. To calculate dqi / d(li / rf) , we, as 
previously, write q\ = Ni/n as (N/N) x (N/n) ~ 
(k/L^Rel = (/ 4 /?7) 4 (ry/ J L) 4 Re| and take into ac- 
count that n/L ~ Re L 3 ^ 4 , which gives Mi(li/n) ~ 
Mi(qi)(k/n) 3 . Also, since the dependence Mj = Mj(gi) 
is not unique (Fig. [1]), the distribution for each given 
Reynolds number Re is taken as the average of the distri- 
butions for Re^ in the range 0.9 < Re^/Re^ < 1.1 with 
the weights of the distributions according to Eq. ([1]). The 
distributions are close to the Gaussian distributions (the 
dashed line curves), with the positions of their maxima 
(in units of the dissipation scale 77) and standard devi- 
ations depending on the Rei, very slightly [as ~ Re^' 07 
and ~ exp(2ReL/10 5 ), respectively]. These results are in 
general agreement with the results of the direct numer- 
ical simulation of isotropic homogeneous turbulence by 
Jimenez et al. [181 ] , who found that the vortex radius dis- 
tributions are bell-shaped, practically do not shift with 
Rei, (with the radius r measured in n units), and have 
maxima at r k 3?y. 

We also note that since the dependence v = v(r) is 
assumed to be the same for all localized structures, the 
velocity space correlation (structure) functions and en- 
ergy spectrum are independent of the size distributions 
of localized structures. With Eq. ([SJ (i.e. with v ~ r 1 / 3 ), 
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they are the same as in the Kolmogorov-Obukhov the- 
ory 0, H, [l9| | . The structure function of p-th order is 
S p (l) — (Svf) ~ Z p / 3 , where I is the space increment, and 
the energy spectrum is E(k) ~ fc -5 / 3 , where k = 2n/l is 
the wave number (see also [3|, l9|). 

The above consideration shows that the proposed 
model captures essential features of turbulent flows. 
Then, an important question to ask is: Why do two 
apparently different approaches, i.e. the flow stability 
analysis with the hydrodynamic (Navier-Stokes) equa- 
tions and the estimate of the statistical weights of the 
collections of localized structures (eddies), lead to sim- 
ilar results? Or, in other words, what is common to 
these approaches to lead to similar results? The statis- 
tical model is based on three assumptions: i) the par- 
ticles are identical, ii) there exists an elementary cell 
of the phase space in which the states of the particles 
are uncertain, and iii) the kinetic energy of the fluctu- 
ations is as in the Kolmogorov theory for the inertial 
interval of scales. The latter assumption seems to be 
not very essential for our purpose. For example, even if 
a quasi-solid rotation of the localized structures is sup- 
posed (in which case Ei ~ N^ 3 ), (q) rapidly increases 
with Re^ as previously, i.e. the statistical preference of 
the structured flow is preserved. The difference is that 
the transition region in this case is located within unrea- 
sonably low Re numbers (4 < Re^ < 50). This is because 
the kinetic energy of the eddies is highly overestimated 
due to the absence of the decrease of the velocity out- 
side the eddy core, which is characteristic of more real- 
istic models, e.g. the Rankine and Lamb-Oseen vortices 
[20( | . The attention should thus be focused on the for- 
mer two assumptions. As well known, the Navier-Stokes 
equation for the gas fluid can be derived from the Boltz- 
mann equation (the Chapman-Enskog expansion [HI), 
which, in turn, can be obtained from the Liouville equa- 
tion (see, e.g., Ref. 13|, l2l|). In the transition from 
the Liouville equation to the Boltzmann equation, first, 
the many-particle distribution function is reduced to the 
single-particle one under the assumption that the parti- 
cles are identical (the Bogoliubov-Born-Green-Kirkwood- 
Yvon hierarchy), and, secondly, because the gas density 
is taken to be low, the interaction between the particles 
is assumed to be a random two-body collision, so that 
the positions and velocities of the particles are uncertain 
within the mean free path A and the thermal velocity 
c, respectively. It follows that both the assumptions are 
implicitly present in the Navier-Stokes equation. A sim- 
ilar rigorous derivation of the Navier-Stokes equation for 
a liquid fluid is problematic because of the many- body 



present model. Since both the approaches, the estimate 
of the statistical weights of the collections of localized 
structures and the flow stability analysis with the hydro- 
dynamic equations, lead to generally similar results, this 
suggests that the driving force for the turbulent transi- 
tion is basically the same as in the present model, i.e. the 
statistical dominance of the structured turbulent state of 
the system over the unstructured laminar state plays a 
key role. The role of the instability is then reduced to 
performing a structural rearrangement of the flow to find 
this state. This statistical aspect of the turbulent tran- 
sition may be essential for a more comprehensive under- 
standing of the nature of this phenomenon. 



mechanism of interaction between the particles [13l. l2l|. 
however, the necessity of the assumptions that the par- 
ticle are identical and the system behaves stochastically 
on the microscale level seems to be evident. 

Resuming, the hydrodynamic (Navier-Stokes) equa- 
tions essentially involve the same two assumptions as the 
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